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We consider the process of primordial black hole (PBH) formation originated from primordial 
curvature perturbations produced during waterfall transition (with tachyonic instability), at the 
end of hybrid inflation. It is known that in such inflation models, rather large values of curvature 
. . . perturbation amplitudes can be reached, which can potentially cause a significant PBH production 

04 ■ in the early Universe. The probability distributions of density perturbation amplitudes in this case 

I can be strongly non-Gaussian, which requires a special treatment. We calculated PBH abundances 

, and PBH mass spectra for the model, and analyzed their dependence on model parameters. We 

obtained the constraints on the parameters of the inflationary potential, using the available limits 
I ■ on Ppbh- 

a; 

^ . PACS numbers: 98.80.-k, 04.70.-s 

(N- 

I \ 

q'- I. INTRODUCTION 

\ According to the observational data (see, e.g., [ij), the primordial curvature perturbation C, is Gaussian with an 

. almost scale-independent power spectrum. It means that the structure of the Universe originated from near-scale 
invariant and almost Gaussian fluctuations. As is well known, in models of slow-roll inflation with one scalar field the 

Q . curvature perturbation originates from the vacuum fluctuations during inflationary expansion, and these fluctuations 

$-H ■ lead just to practically Gaussian classical perturbations with an almost flat power spectrum near the time of horizon 

jj^ , exit, in a full agreement with the data. So far, there is a weak indication of primordial non-Gaussianity (at (2 — 3)a 

■ level) from the Cosmic Microwave Background (CMB) temperature information from the WMAP 3-, 5-, 7-year data 

^: ii- 

■ It has beenpointed out long ago that for inflation with multiple scalar flelds possibilities exist for non-Gaussian 
I> I fluctuations In particular, authors of @ had elaborated a model of Cold Dark Matter (CDM) (motivated 

. by double inflation scenarios) in which it had been assumed that the initial perturbation field (the gauge-invariant 
" potential) is a combination of a Gaussian field 4>i and the square of an another Gaussian field 02 , ^' = (Z*! + (^2 1 and, 
besides, that (j)2 is described by a sharply peaked power spectrum. 

The detectable non-Gaussianity is predicted in models with additional scalar flelds contributing to C- The time 
CN I evolution of the curvature perturbation on superhorizon scales (which is allowed in double inflation [7] and, in general, 
] ■ in multiple-fleld scenarios) implies that, in principle, a rather large non-Gaussian signal can be generated during 
^ I ' inflation. The primordial non-Gaussianity of C in multiple-fleld models can be calculated using the SN approach 
or the expression for ^ through the non-adiabatic pressure perturbation [TT|-[l3|. It is important to note that 
^ ' non-Gaussian contributions to C, predicted by all these approaches might be negligible on cosmological scales but 
k>( , rather large on smaller scales allowing, in principle, primordial black hole (PBH) formation. 

■ There are several types of two-field inflation scenarios in which detectable non-Gaussianity of the curvature pertur- 
d [ bation C can be generated: curvaton models [l3 - [T9j , models with a non-inflaton fleld causing inhomogeneous reheating 

[13, HH, curvaton- type models of preheating (see, e.g., j22| and references therein), models of waterfall transition that 
ends the hybrid inflation [23l-[33j. 

In these two- fleld models, the primordial curvature perturbation has two components: (^g, which is a contribution 
of the inflaton (almost Gaussian) and Co-, which is a contribution of the extra fleld a. This second component is 
parameterized by the following way (33 |: 

C.(x)=aa(x)+a2(x)-(a2). (1) 

If the linear term in ((IJ is negligible (i.e., if a « 0), one has the "x^-model", in which curvature fluctuations are 
described by ^^-distribution. This x^-model is a particular case of Xm-niodel (35l - [37| . i.e., a model in which the 
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fluctuations generated during inflation have x^,j-distributions. Of course, there is a severe observational constraint 
on the spectrum amplitude V,^^ [38] predicted by the x^-model, at cosmological scales. That is, the dominance of 
the quadratic term in ([T]) and, correspondingly, large non-Gaussianity are possible only on smaller scales (where, in 
particular, in case of a blue tilt the amplitude Vq^ can be of order of one, leading to PBH formation). The possibilities 
of PBH formation in curvaton-type scenarios are discussed in (s^ . 

It was shown in the previous work of authors [s^ that the power spectrum of the primordial curvature perturbations 
from the waterfall field in hybrid inflation with tachyonic preheating has a form of a broad peak, and the peak value, 
fc*, depends on the parameters of the infiationary potential, in particular, on the parameter /?, which is the ratio 
|m^|/iJ^, where is the mass-squared of the waterfall field x- At small values of /3 (/3 < 10) the peak is far beyond 
horizon (k^/aH ^ 1) and the perturbations are strongly non-Gaussian (because they have ^^-distributions due to 
the fact that curvature perturbation C depends on the waterfall field amplitude quadratically, just like in Eq. ([1])). It 
appears that the spectrum amplitude is negligible on cosmological scales but is quite substantial at small scales 
(if /3 ~ 1) and it is interesting to analyze if it can be constrained by data of PBH searches. 

The effects of non-Gaussian primordial curvature and density perturbations on the formation of PBHs had been 
considered in works (40l - |4^ . Our consideration in the present paper has almost no intersections with the content 
of these works. We study in this paper, mostly, two questions: i) forms of the PBH mass spectra in non-Gaussian 
case and ii) constraints on the inflationary potential parameters (for the concrete inflation model) following from 
processes of PBH formation in radiation dominated era. Note that we consider the case of the strong non-Gaussianity 
(in contrast with, e.g., studies of 44]). 

The plan of the paper is as follows. In Sec. [Illwe review, very briefly, main aspects of the model used for a describing 
of the waterfall transition (the tachyonic preheating) at the end of hybrid inflation. In particular, we present in this 
Section the formula for the curvature perturbation ( (derived in our previous work [33|), which is a basis for all 
following calculations of PBH formation. In Sec. IIIII we consider problems connected with the process of PBH 
formation in radiation era of the early Universe: probability distributions functions of our model, formula for PBH 
mass spectrum (following from Press-Schechter formalism), formula for the relative energy density of the Universe 
contained in PBHs. At the end of Sec. IIIII we present some illustrative results of PBH mass spectra calculations. In 
Sec. IIVI we present the resulting constraints on the parameters of our inflation model following from the studies of 
PBH formation. Sec. IVl contains our conclusions. 



II. THE WATERFALL TRANSITION MODEL 



We consider the hybrid inflation model which describes an evolution of the slowly rolling inflaton field (p and the 
waterfall field x, with the potential (45|, |46| 

^(0, X) = - ^X'^ + lrr^'cb' + (2) 

The first term in Eq. is a potential for the waterfall field x with the false vacuum at x = and true vacuum at 
Xo = 2M^/\/X = v"^. The effective mass of the waterfall field in the false vacuum state is given by 

ml{ct>) = J c^l) , (3) 

At (jy^ > (jy^ the false vacuum is stable, while at 0^ < (/)^ the effective mass-squared of x becomes negative, and there 
is a tachyonic instability leading to a rapid growth of x-modes and eventually to an end of the inflationary expansion. 
The evolution equations for the flelds are given by 

X + Si/x-V^x = (2\/AAf 2 - _ Ax2)x. (5) 
From Eq. ([5]), one obtains the equation for Fourier modes of 6x'- 

dxk + 3HSxk + - I3H! + T'/-') Sxk = 0. (6) 
Here, a = a{t) is the scale factor and the parameter /3 is given by the relation 

/\,f2 

/3 = 2VA^. (7) 



3 



The solution of Eq. ^ (in which we ignore gradient term due to the choice of a uniform 0-gauge) is (for t > tc, tc 
is the critical point when the tachyonic instability begins) 



The time evolution of Sxk during the growth era of the waterfall was studied, numerically, in the previous work of 
authors [1^. We used in [SJI an artificial cut-off of large-fc modes, which corresponds to considering only the waterfall 
field modes that already became classical near the beginning of the growth era. The classical nature of the waterfall 
field had been discussed in (30j . in the approximation when the expansion of the Universe is negligibly small. It has 
been shown in 30] that, in the Heisenberg picture of the quantum theory, the operator Sxk has, at not very large fc, 
almost trivial time dependence (during the most part of the growth era) . Namely, Sxk is a constant operator times 
a c-number, which means that the perturbation is classical (this issue is elaborated in detail in the literature on the 
quantum-to-classical transition [47h49| ). 

Following [s^] we assume that the waterfall transition ends when the last term in right-hand side of Eq. (|4]) becomes 
equal to the preceding one, i.e, when 

{iSx?)^ — ^xli- (9) 
7 

The main equation for a calculation of the primordial curvature perturbation (on uniform density hypersurfaces) 
is [TH (see also fr2, 13]) 

C = - / dt^^^^ (10) 
J P + P 

where the non-adiabatic pressure perturbation is Spnad = ^E,~ c^s^P ^^.d the adiabatic sound speed is = p/ p. The 
formula (|10p follows from the "separated universes" picture [7Hll| where, after smoothing over sufficiently large scales, 
the universe becomes similar to an unperturbed FRW cosmology. In our case, one has 

P 

Spnad ^ 5p^ - 5px- (11) 

Energy density p and pressure p is a sum of contributions of 4> and x fields. Eq. ([Tl]) takes into account that in 5pnad 
there is no contribution from (j) field. 

For the curvature perturbation, we have the formula 

C = Cx = -^(x'-(x')) (12) 
where x^ and (x^) are determined at the time of an end of the waterfall, t — tend, and A is given by the integral 



tend . ,2 



A = 



Hcdt f fit) Y 1 



<P'{t) + {x'{t)) \fitend) 



fit) / P 




(13) 



Here, the function f{t) describes the time evolution of the waterfall field, which is almost independent on k |30l. |33|. 

x(x,t) = C(x)/(t). (14) 

It was shown in [sl] that for /3 ^ 1, the curvature perturbation spectrum will reach values of ^ 1 in a broad 
interval of other model parameters (such as r, 7, He). The peak values, fc*, for small f3, are far beyond horizon, 
so, the smoothing over the horizon size will not decrease the peak values of the smoothed spectrum. Furthermore, 
the spectrum near peak remains strongly non-Gaussian after the smoothing. The calculations of [33| . based on the 
quadratic inflaton potential, show that for /? ^ 100 and in the broad interval of r the peak value k^ can be estimated 
by the simple relation: 

where N is the number of e-folds during the waterfall transition. The similar estimate is contained in the recent work 
[3^ 1. Note that for k < fc*, we obtain the well-known (e.g., [Ullill^l) result: 'P^(fc) ~ k^. 
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III. PBH PRODUCTION FROM NON-GAUSSIAN PERTURBATION 
A. PBH formation threshold 



A production of PBHs (about these objects, see, e.g., reviews [50ll5l|) during reheating process had been considered 
in works 52 56]. PBH formation in connection to non-Gaussianity has been studied in 40-44i]. In ah those papers, 
the case of rather weak non-Gaussianity has been considered. In the present work, we study in detail the case of 
strong non-Gaussianity (i.e., one when quadratic term in Eq. ([I} dominates). Furthermore, we have an opposite sign 
for the quadratic term due to sign in Eq. (fT2|). which leads to very different dependencies of PBH abundances on the 
power spectrum amplitude cornpared to the usually considered (Gaussian or almost Gaussian) cases (such behavior 
was qualitatively described in [30|). 

The classical PBH formation criterion in the radiation-dominated epoch is [H^l 

5 > (5c w 1/3, (16) 

where S is the smoothed density contrast at horizon crossing (at this moment, k — aH). The Fourier component of 
the comoving density perturbation 5 is related to the Fourier component of the Bardeen potential as 



2 / fc ^ ^ 



(17) 



For modes in a super-horizon regime, ~ — (2/3)Cfc, so (|16|) can be translated to a limiting value of the curvature 
perturbation [isj . which is 

Cc = ^ 5c « 0.75. (18) 
If we assume somewhat larger PBH formation threshold, 5c ~ 0.45 (see, e.g., [1^), then 

Cc « 1. (19) 

We will not insist on the concrete value of the threshold parameter and, in the following, will consider both values 
([T8| and (fT9|) as possible ones. 

B. Perturbation probability distributions 

The relation between curvature perturbation C and the waterfall field value is given by Eq. (|12p. or, using cr^ — (x^), 

C = -Aix^ - = Crnax - Ax\ (ma. ^^a^. (20) 

Here, A and cr^ generally depend on the smoothing scale R. The distribution of x is assumed to be Gaussian, i.e., 

2 

Pxix) = ^^e'^. (21) 



The distribution of ( can be easily obtained from (^01 [^T|) : 



dx 
dc 



eVjir, C<Cmax, (22) 



which is just a x^-distribution with one degree of freedom, with an opposite sign of the argument, shifted to a value 
of Cmax- As required, {() = and 



(C')= / C'pc(C)rfC = 2CL.- (23) 



On the other hand. 



(24) 
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FIG. 1: The probability distribution Pq{C) for E = 0.7, P° = 0.4 and for different values of 7?: from left to right, R — 
10kg ^ , kg^, 0.1kg ^ . TliB posslMe (model dependent) values of (^c (Cc = 0.75 and Cc = 1) are shown by dashed lines. 

where W(kR) is the Fourier transform of the window function, and we use a Gaussian one, W^{kR) = exp(— fc^i?^), 
in this work. 

From 1221 1211) we can write for (max (we now denote the argument R explicitly): 



Cmax {-^) 



I j Vdk)WHkR)f 



1/2 



Everywhere below, we use the following notation: Qmax{R — 0) = Qynax- So, 



1 f ^ ,, .dk 



1/2 



(25) 



(26) 



It is clear that PBHs can be produced in the early Universe, if C,max > Cc 



C. Formulae for PBH mass spectrum and abundance 

In general, a Press-Schechter approach [s^ used for calculations of the number density of clusters in different 
scenarios for the formation of structure in the Universe, does not assume a use of the assumption that the initial 
perturbations have just Gaussian distributions. For example, in [60j the Gaussianity of the cosmological density field 
was tested using two different models for probability distribution functions: a standard Gaussian model and a texture 
(see, e.g., [6l|) model. In [l^ the Press-Schechter formalism had been used for the case when an initial density field 
has a x^-distribution, and PBH abundances, as a function of black hole mass, for a power-law primordial power 
spectrum, were calculated. 

The energy density fraction of the Universe contained in collapsed objects of initial mass larger than M in Press- 
Schechter formalism [s^ is given by 



1 

Pi 



Mn{M)dM = / pc(C)dC = ^'(C > Cc;i?(M),i^), 



(27) 



M 



Cc 



where function P in right-hand side is the probability that in the region of comoving size R the smoothed value of C 
will be larger than the PBH formation threshold value, n{M) is the mass spectrum of the collapsed objects, and pi 
is the initial energy density. Here we ignore the dependence of the curvature perturbation C, on time after the end of 
the waterfall, assuming it does not change in super-horizon regime, until the perturbations enter horizon at A: = aH . 
The mass spectrum of the collapsed objects, n(M), is given by 



n(M) = 2^ 



dP 
dR 



dR 



(28) 
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where, as usual, the factor 2 approximately takes into account the fact that underdense regions also collapse. In the 
above formula, M is the initial fluctuation mass corresponding to the fluctuation of the comoving scale R, 

M = ^Ma^Rf; ^ = (4.)-V33-2/3^-i/3„-i^,^-2/3 (29) 

(fli is the value of the scale factor at t ~ ti, M is calculated at the moment ti corresponding to the time of the end of 
the waterfall, which is assumed to be close to reheating time; note that comoving fluctuation mass is not constant). 
The horizon mass corresponding to the time when fluctuation with initial mass M crosses horizon is (see 62]) 

Mh = M^/^m2/3, (30) 

where Mj is the horizon mass at the moment ti, 

47r , An, ,,0 AttMI , , 

« —t\p^ « y (H-^)V = (31) 

(here, we used Friedman equation, pi = 3MpH^). The reheating temperature of the Universe is [62j 

/90M|i/^y/* 



Trh = ( ^/^ ' ) , « 100. (32) 



For simplicity, we will use the approximation that mass of the produced black hole is proportional to horizon mass, 
namely, 

Mbh - fhMh - /,,Af//'M2/3, (33) 

where fh ~ (1/3)^/^ — const (this particular value of fh corresponds to a threshold of the PBH production in Carr- 
Hawking collapse, see, e.g.. Appendix of paper ^^). Our final qualitative conclusions do not depend on the value 
of fh- In more accurate analysis, one must take into account that the connection between Mbh and Mh is more 
complicated, and, in particular, depends on the type of the gravitational collapse [13, [H, El, [gJ]. Moreover, there 
can be a dependence on a shape of the radial fluctuation profile [63 |. 

Using Eqs. (|281 [29l [33|) . the PBH number density (mass spectrum) can be written as 



nBH{MBH)=n{M) 



dM /47r\~^/^ dP fhp]'^Ml'^ 



cLMbh \ 3 



dR 



aiMljj 



(34) 



We can estimate the relative energy density of the Universe contained in PBHs, at the moment of time t (assuming 
radiation-dominated Universe with t < t^q and ignoring the PBH mass change due to accretion or evaporation) as 
follows 



1 / \ 3 
l a 



^PBH{t) W Xaitj) J ^^BHMBHdMBH- (35) 

Using the scaling relations p ^ a ^ i^/^ and considering the moment of time for which horizon mass is equal to 
Mh, we obtain 

1 / M \ ^^"^ f 

VLpBH{Mh) ~ - i-j^j j nBEMl^dhiMBH. (36) 

It is well known that for an almost monochromatic PBH mass spectrum, D,pBH{Mh) coincides with the tradition- 
ally used parameter (3pbh (energy density fraction of the Universe contained in PBHs at the moment of their for- 
mation). Although all PBHs do not form at the same moment of time, it is convenient to use the combination 

Afj ^^"^ p~^ Mg^^nBHiMBn) to have a feeling of how many PBHs actually form. We will use this combination in the 
following (Fig. |3] below). 

One should note that, strictly speaking, observational constraints on the parameter /3pbh are obtained using the 
assumption that the PBH formation takes place at a single epoch and the corresponding spectrum has a narrow 
peak at some value of k. In our case, a width of the peak is not very small (see Fig. 4 in 33]). But, as is well known, 
the black hole abundance is extremely sensitive to an amplitude of P^(fc) and, correspondingly, the bound on 'P^(fc) 
weakly depends on the value of Ppbh- Our main aim is to determine constraints on the inflation model parameters 
(in particular, on the value of /?), and, as we show in the present paper, the black hole abundances depend on /3 
and on amplitude of V^^k) very strongly. In such a situation the form of the T^i^-spectrum in real scenario of PBH 
formation (a size of the spectrum's width) is not very essential. 
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FIG. 2: The probability P{C > (c) as a function of 7?. From top to bottom, P° = OAiCmax ~ 0.9), 0.28(0 
0.752), 0.27846(Cma:i ~ 0.750012). For all curves, S = 0.7, Cc = 0.75. 



D. PBH mass spectrum calculation 

The power spectra V(^{k) from waterfall transition process, for different sets of parameters, have been calculated in 
our previous paper 33] . For the purpose of this section, it is convenient to parameterize the curvature perturbation 
power spectrum as follows, 



■p^(fc) =750 exp 



(Igfc/fco 



,21 



2S2 



(37) 



where gives the maximum value approached by the spectrum, fco is the comoving wave number corresponding to 
the position of the maximum and S determines the width of the spectrum. Evidently, values of the parameter ka 



should be equal to the corresponding peak values, A;*, of the "P^ (/c)-curves calculated in [33|. 

What are the typical values of the parameters and E? It was shown in [s^l that the spectrum V(^{k) can reach 
values of order of 1 if /3 1. It follows from Fig. 4 of [H] that for /3 = 2, r = 0.1 one has E w 0.7, « 0.4 (and 
Cmax ~ 0.9), while for /3 = l,r = 0.1 the peak is more wide and high: E « 1.0, « 1.2 (with (max ~ 1.86). 

So, in case of /3 = 2 the PBHs are produced if formation threshold is assumed as in ((T5)l while they are not yet 
produced in case of (|19p. 

We note that, generally, the values of V^, E and fco depend on the rather complex interplay between the potential 
parameters /3, r, Hj and 7. In particular, the position of the peak, fco, depends on Hj and the number of e-folds 
which waterfall takes, and thus varies with the change of any of the above 4 parameters. So, in this Section, we will 
consider 7^°, fco and E as independent parameters. 

In Fig. [T]we show the distribution P({C) (Eq. ([H])) for a fixed set of parameters E, (physically, they correspond 
to /3 — 2, r — 0.1 case) and for different values of the fluctuation size (smoothing scale) R. It is seen that in this case, 
for larger values of R the fluctuation spectrum is not "strong" enough to produce PBHs, while for smallest R the 
PBHs will form (if we assume threshold ((T5)) - both values of Cc that we consider are also shown in the Figure). 

The probability P{( > Cc) (used in Eq. (P7)) ) for fixed values of E and Cc, but for different is shown in Fig. 
The value of was fine tuned for the bottom curve so that (max does not go below Cc (and (max — Cc ^ !)■ Once 
Cmax drops below Cc no PBHs ever form (at least for the classical PBH production scenario which we consider). For 
the bottom curve, it happens at fco-R <C 1 just because Cmax{R) — Cc > only for the smallest values of R (see Eq. 
(p5)) ). For each curve in Fig. [2]we see a sharp drop of P(C > Cc) to the zero at the value of R when (max{R) reaches 
Cc. Technically, the derivative dP/dR at this point diverges, which will lead to a characteristic spike in the PBH mass 
spectrum, according to Eq. ([M)) . 

Results of the PBH mass spectra calculations using formula (p4l) are shown in Fig. [3] Again, we are doing some 
fine-tuning for part of the curves, choosing parameters so that (max — Cc ^ 1- In these cases, such fine tuning allows 
us to reach Ppbh ^ ^^'^ p^^ Mg^nBH{MBH) ^ 10"'^ or so. The mass of the produced PBHs, as can be seen from 
the same Figure, is several orders below M^. 
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FIG. 3: The PBH mass spectra for different values of the perturbation spectrum amplitudes. From right to left, Vq = 
l{(,ma:c ~ 1 .42) , 0.4(C™a.= « 0.9) , 0.28(Cma. « 0.752) , 0.27846(Cma. ~ 0.750012) , 0. 2784511 (Cma. ~ 0.750000068). The position 
of the peak in P^(fe)-spectrum is the same for all cases. For the calculation we used the value E = 0.7, and C,c = 0.75. The mass 
corresponds to horizon mass at the moment of time when perturbation with comoving wave number fco enters horizon. 

It is well known 51] that for rather large PBH masses, say, Mbh ^ 10^" g, the constraints on their abundance 
from different types of sources (such as gravitational constraints and constraints following from non-observation of the 
products of PBH's Hawking evaporation) are rather severe, (3pbh ^ (10~^^ 10"^'^), or so. It turns out that in the 
present model it is possible to reach such small values of Ppbh only with extreme fine tuning of inflationary potential 
model parameters. However, in the range of masses below Mbh ~ 10^" g the constraints on Ppbh are not so severe, 
and in this mass range, Ppbh ~ 10"'^ is not forbidden by the observations. Such light PBHs evaporate very lastly 
(long before nucleosynthesis) and the hope of their possible detection is mainly due to high-frequency GWB which 
they generate through Hawking evaporation (see, e.g., [65l - l67j : in these papers, PBHs with masses Mbh ^ 10^ g and 
values of fipBH ~ 10^"^ are considered). 

IV. CONSTRAINTS ON THE WATERFALL TRANSITION MODEL 

We have seen that in the waterfall model considered, PBH abundance severely depends on the amplitude of the 
curvature perturbation spectrum: once Cmax is above Cc, PBHs are produced intensively. Demanding that PBHs do 
not form in the early Universe, we can impose the bound on parameters of the inflaton potential. From the condition 
Qmax < Cc one has, for two fixed values of Cc, the following constraints (r = 0.1; there is a weak dependence on this 
parameter but the result is almost independent on j,Hc)'. 

Cc = 0.75 : /3 > 2.3, < 0.29; (38) 

Cc = 1 : 13 > 1.65, < 0.55. (39) 

These constraints are not based on the comparison with data on Ppbh- 

It is interesting to estimate the mass of PBHs that can be produced by this model, and corresponding horizon 
masses. For the horizon mass corresponding to the peak position ^33)] one has the relation 

Ml « e^^Mi, (40) 

where N is the number of e-folds that waterfall transition takes. This approximate relation follows from the estimate 
k^/aH ^ obtained in [33| . Here, fc* is the peak value of the C power spectrum, aH corresponds to an end of the 
waterfall. We find that N is highly parameter-dependent: although mostly depends on /3 and rather weakly on 
other parameters, N depends on 7 (or 0c) in a strong way. Here we estimate the range of N values to get an idea 
what PBH mass range is given by Eq. PO)) . 

For the purpose of this estimate, we take the maximal possible value of (j)c to be equal to the (reduced) Planck 
mass Mp. The lower limit is obtained from the classicality condition [30]: to have an effective waterfall (classical 
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FIG. 4: The horizon mass regions corresponding to the position of the peak in curvature perturbation power spectrum (shaded 
areas), a) /3 ~ 2.3, ~ 0.75; b) /? ~ 1.65, = 1- In both cases, the value of P is just enough to produce PBHs, and we 
vary the parameter <l)c (or, equivalently, 7) between minimal and maximal possible values. The exact PBH abundance and 
relation of characteristic PBH mass Mbh to Mh will depend on values of parameters in a fine-tuning regime (see Fig. |3]for an 
illustration) . 



regime is reached for the x field), we need ^/J ^ ^/VP, or, at least, 7 ~ 1//3 (numerical solutions show that for sets 
of parameters we consider, the waterfall is still effective in this case). So we consider the lower limit for (j)c to be 

We find that w 4 for the limit of 0c = /3-ffc, both for f3 = 1.65 and /3 = 2.3; this result turns out to be independent 
on He- With the growth of 0c, N also increases. For 0c = Mp and He = 10" GeV, TV w 32 for /3 = 2.3 and iV sa 41 
for /3 = 1.65 (these values of N also have some dependence on inflation energy scale He). 

We plot the possible regions of Mh in Fig. |4l It is seen that the possible mass range is very wide. As discussed 
above, the exact PBH abundance (Ppbh) and relation of characteristic PBH mass Mbh to Mfi will depend on values 
of parameters in a fine-tuning regime (see Fig. [3] for an illustration). 



V. CONCLUSIONS 



We have considered PBH production from strongly non-Gaussian density (curvature) perturbations in the radiation- 
dominated era of the early Universe. The main physical model that we used is the model of hybrid infiation waterfall 
(with tachyonic preheating), however, the results may be applied to different models that produce similar perturba- 
tions. 

We have given the expressions for perturbation probability distributions and, on the basis of Press-Schechter 
formalism, calculated PBH abundances and PBH mass spectra for the model. The most important result of the paper 
is that we obtained limits on the parameters of the potential of our hybrid inflation model. In particular, we have 
shown that the parameter /3, which is the ratio jm^l/ff^, is limited from below, i.e., /3 is larger than some value 
(otherwise the abundance of PBHs will be too high, in contradiction with available constraints on Ppbh)- We have 
shown also that these limits on /3 are sensitive to the PBH formation threshold parameter (in our case, Cc)- Note that 
the characteristic PBH masses that can be (in principle) produced by this model, are shown to depend signiflcantly on 
the coupling parameter, 7, of the inflaton potential. The second important result is that for our inflation model the 
possible horizon and PBH mass regions, corresponding to the peak in perturbation power spectrum, are constrained 
(Fig. H. 

It was shown also, that to obtain values of PBH abundance not contradicting with the available limits on Ppbh, 
for Mbh ^ 10"'^° g, extreme fine-tuning of the model parameters are needed. Note that the model allows to more or 
less naturally (with a degree of fine-tuning similar to one needed in some single- field infiation models [44, 68]) produce 
PBHs with Mbh ^ 10^ g or so, and (3pbh ^ 10^"^ or so. Such PBHs are one of the possible sources of high-frequency 
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GW background, which can be observed in future. 

Throughout the work, we used simphfied gravitational collapse model (using Eq. ([55]) ). close to the standard one, to 
treat PBH formation. The inclusion of critical collapse [H, l63l | effects in our calculation would require the replacement 
of fii in Eq. ([55]) with a function proportional to — Cc)^"! 7c is around 0.3 — 0.4. This will change the results for 
PBH mass spectrum (Fig. [3]) - in particular, the curves on that Figure will be shifted by 1 — 2 orders of magnitude 
to the left and their shape will change, but this will not affect other results of the paper. 

Note added. After this work has been published as an e-print, the paper [69| appeared, in which an analytical 
framework for calculating the curvature perturbation spectra produced by the hybrid inflation waterfall is developed, 
in a general case of any inflaton potential form including the case of N > 1. 

The author of (69j finds good agreement between numerical calculations of our work (ssj and his analytic estimates. 
He also discusses in detail all assumptions that are made in such models, in particular, he stresses that the gradient 
of Y must be negligible for the calculation using formula PH)) to be viable. We proved the smallness of this gradient 
in [sl] for our case (see Eq. (3.9) of that work). Also, discusses the compatibility of the evolution equation for x 
(our Eq. ^) with the energy continuity equation. The possible inconsistency between these two equations is due to 
the fact that interaction term, ^ 0X^i is dropped in Eq. (H]) but term ~ (jx^x is still used in Eqs. ([5l[6|). The condition 
for using this approximation correctly is derived to be [69| 



In our case, from Eqs. Q and ([S]), one has 



I ^!%M«i. (41) 



(42) 



which is small only for A'' = Ht > 1. For smaller Ht, when the waterfall just started, Xk is small, and term proportional 
to (px^ in Eq. (H]) can be dropped just due to smallness of x- As one can see from Eq. (j35]), the consistency condition 
depends only on a value of the parameter r which we fixed throughout the present work. Namely, we used the value 
r — 0.1 and, in this case, the consistency condition is satisfied. 
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